Afc# at- 


NASA-CR- 177362 
19850022732 


A Reproduced Copy 

OF 


Reproduced for NASA 
by the 

NASA Scientific and Technical Information Facility 


LIBRARY COPY 


\\U / 1 


1:485 


LANGLEY RESEARCH CENTER 
LIBRARY, NASA 
HAV.nON, VIRGINIA 


FNo 672 Aug 65 


NF00675 



1 Report No J 2. Go*irr*n®r>i No 

NASA CR- 177362 | 

3 R*Cip*tnl | CeUKrJ hLx 

4 Titu W'd SuMlOa 

Estimation of Dynamic Rotor Loads for the Rotor 
Systems Researcn Aircraft: Methodology 
Development and Validation 

5- Report D*tt 

Mav, 1985 

6. Fir#orm>r*j G«^rw*t>o*"i Qxf* 

FHF 

7 Awffortl) 

R. W. Du Val and M. Bahrami 

8 Perform k>j CVpieutttvOA Report No 

13 Work Umt No. 

Job Order # T4713D 

4 Pwforrr»*$ Adtfreti 

Advanced Rotorcraft Technology, Inc. 
1804 Stierlin Road, Suite 210 
Mountain View, CA 94043 

11 Co^i/krt or Oent 

NAS2-1 1688 

13 Type of Report KVj r«f«od Cov*r »0 

Contractor Report 

12 Spomonr*} AgstYKY Aqdmi 

National Aeronautics and Space Administration 
Ames Research Center, Moffett Field, CA 94035 

14 A#** cy Gx* 

RT0P # 505 4251 


IS Sw0O**«vmnrv Nom 


Point of Contact: Technical Monitor, C. W. Acree, M.S. 237-3, 

NASA Ames Research Center, Moffett Field, CA 94035 

(415) 694-5664 or PTS 464-5664 ; 

16 AtrtTACt 

The Rotor Systems Research Aircraft uses load cells to isolate the 
rotor/transmission system from the fuselage. A mathematical model relating 
applied rotor loads and inertial loads of the rotor/transmission system to 
the load cell response is required to allow the load cells to be used to 
estimate rotor loads frcm flight data. In this study, such a model is 
derived analytically by applying a force-and-.mo.rent balance to the isolated 
rotor/transmission system. The model is tested by comparing its estimated 
values of applied rotor loads with measured values obtained from a ground 
based shake test. Discrepancies in the comparison are used to isolate 
sources of unmodeled external loads. 

Once the structure of the mathematical model has been validated by 
comparison with experimental data, the parameters must be identified. 

Since the parameters may vary with flight condition it is desirable to 
identify the parameters directly frcm the flight data. A Maximum Likelihood 
identification algorithm is derived for this purpose and tested using a 
computer simulation of load cell data. The identification is found to 
converge within 10 samples. The rapid convergence facilitates tracking 
of time varying parameters of-the load cell model in flight. 


17 Kty Word* iSu&tnri fcy Aut*or(ilI 

Aeronautics, Load Cell, Dynamic, 
Calibration, Rotor Systems Research 
Aircraft, Parameter Identification 

13 0‘«rit*.~icn 

Unclassified - v limited 
Subject Category 05 


10 Srctxity Oreif (of l*vt report) 

Unclassified 

20 SacvrtTY CXrtfcf (of ttwi F ttyzi 

Unclassified 

| 21 cf 

50 

22. 


*For ut* t*i# Nation*! Ttc*ifucil C*rv»ci Vir^m* 22101 



NASA CONTRACTOR REPORT 


177362 


JHASA-CH- J77362) ESTIMATION OF DYNAMIC 
aoroa LOADS FOB THE BCTOB SYSTEMS BESEABCB 
AIRCRAFT: HETDCDCLGG Y DEVElCf BENT AND 

VALIDATION j^Advauced fotorcraf t Technology, 

Inc.) 5 ft F He A04/KP AOt CSCL 01C G3/05 


H85-3 1045 

UncJaa 

24912 


I 


Estimation of Dynamic Rotor Loads for the Rotor 
Systems Research Aircraft: Methodology 
Development and Validation 


R. W. Du'Val 
M. Bahrani 


CONTRACT NAS2- 11683 
MAY 1965 



y 

AUG1S35 

RtCEiVEO 




NASA CONTRACTOR REPORT 177362 


Estimation of Dynamic Rotor Loads for the Rotor 
Systems Research Aircraft: Methodology 
Development and Validation 


R. W. Du Val and M. Bahramj 
Advanced Rotorcraft Technology, Inc. 
Mountain View, California 
May, 1985 


Prepared for 

Ames Research Center 

under Contract HAS2-116S3 * 


Nancnal Aeronaut, cs and 
Space Acm.nistrat-on 


Ames Research Center 
Mc'fett Fieto Ca | i , c r n>a 94035 



TABLE OF CONTENTS 


Section Page 

Introduction . . . , . . . l 

Nomenclature . . . . , . . , 2 

Phd "0 I Development and Aalidation of an Analytical model of the Rotor 


Tran<rnibion Load cell system ... . a 

Background . . 5 

Preyious Methods . . . . , . . ... 5 

Proposed Method » * *. ^ . ... 5 

Objective and Approach . . . . . . G 

Rotor Transmission Model . ... ... 6 

Test Conditions . . .... 7 

Error Analysis .... . . . 7 

Total Estimation Error . 7 

Systematic Errors In Data Collection . . . 8 

Inmodeied Static And Dynamic Effects . ... 9 

Random Errors In Sensors . . . 9 

Lnrnodeled External Loads . . . . ... . . 10 

Incorrect Parameter A’alues . 10 

Results . . ... 10 

Conclusions For Phase l . . II 

Concluding Pemarks . . 11 

Specific Conclusions . . 12 

Recommendations For Further Re>eurh . 12 


//fs- * r_ 


Phase II Identification of Rotor/Transmisston Load Cell model parameters direct!) 


from flight test data . 

Background . , . , . . 

Overview . ... 

Identification And Calibration . . .... 

Objective And Approach 

Design Of The On-Line Calibration Algorithm 
Approach ... 

The Load Cell Model . ... 

A Dynamic Model of Applied Loads . . 

Design of The Kalman Filler 
The Maximum Likelihood Estimation Algorithm 
Operational Approach . 

Simulation Results 

Identification From Ambient Vibration^ 

Conclusions For Phase II 


14 

u 

14 ' 

14 
16 
16 
16 
17 

15 
19 
21 


24 

24 


Concluding Remarks 
Specific Conclusions 
Recommendations For Further Research 
References . ... 

Appendix A Load Cell Model Derivation 

Appendix B The effect of Sensor Noise on Load Estimation 

Appendix C. Calibration Mgomhm 

Figures . ... 


24 


30 

31 

»>*> 
KJ m 


II 



! 


I a 


LIST OF TABLES AND FIGURES 


Table I * Comparison of Actual Estimation Levels with Predicted Levels Based on Instrumentation 
Noise 

Figure 1 Loads acting on isolated rotor/transmission 

Figure 2 Load cell arrangement for rotor/transmission svstem 

Figure 3 RSHA rotor balance s>v stem and mounted itemitems. 

Figure \ Magnitude of load cel! transfer functions for 800 lb applied lateral force 

Figure 5 Magnitude of inertial load transfer functions for 800 lb applied lateral force 

Figure G Magnitude of x-axis applied force estimation error for applied lateral forces 

Figure 7 Magnitude of y-axis applied force estimation error for applied lateral forces. 

Figure 8 . Magnitude of r-axis applied force estimation error for applied lateral forces 

Figure 9 Magnitude of x-axis applied moment estimation error for applied lateral forces 

Figure 10 Magmtud of y-axjs applied moment estimation error for applied lateral forces 

Figure 11 Magnitude of x- axis applied moment estimation error for applied lateral forces 

Figure 12 Comparison of magnitudes of x-axis force estimation error for applied x. y and z axis 

loads 

Figure 13 Comparison of magnitude of y axis force estimation error for applied r y and z axis loads 

Figure 14 Comparison of magnitude of c-axis force estimation error for applied x y and x axis lo'ads 

Figure li Comparison of magnitude of x-axis moment estimation error for applied x, y and r-axis 
loads. 

Figure 16 Comparison of magnitude of y-a\is moment estimation error for applied x, y and x axis 
loads 

T igure 17 Comparison of magnitude of x-axis moment estimation error applied x, y and x-axis loads 
Figure 18 Phase lag of unmodeled roll moment for applied lateral loads 
Figure 19 Schematic diagram of on-line calibration 
Figure 20 Identification of vertical load cell angle 
Figure 21 Identification of lateral load cel! angle 
Figure* 22 Comparison of actual and estimated y-axis hub loads 
Figure 23 Comparison of actual and estimated rolling moment hub loads 
Figure A 1 Undcformed load cell geometrv with generalized reaction force 
Figure A 2 Load tell geometrv with generalized reaction force components 

m 




? 

?>. 


i 


t 


i 

4 

i 

i 


! 






INTRODUCTION 


The Rotor Systems Research Aircraft (RSRA) has a set of seven load cells connecting the mam rotor 
transmission to the fuselage Their purpose is to make high accuracy measurements of the net rotor loads 
as resolved at the rotor hub from flight data 1 . The use of these load cells to estimate applied rotor 
forces and moments at the hub requires an accurate mathematical expression relating rotor loads, inertial 
loads and load cell readme Both the structure and parameters of this model must be specified Previous 
approaches ut processing ground test data have not yielded an acceptably accurate relationship for the 
ca^o of high frequency dynamic lo«tds 

The previous approaches have utilized empirical techniques ba^ed on experimental data to obtain the 
desired model The approach presented in the first phase of this study to develop the model analytically 
and validate it from experimental data using modern techniques of model structure determination and 
parameter identification This approach provides a level of physical insight that has been missing in 
the purely experimental approaches and the resulting analytical relationship provides a unified basis for 
utilirmg both static and dynamic calibration data to validate and upgrade a common model 

The second phase of this study develops and demonstrates a technique for identifying the parameter 
values of the model directly from flight data thereby eliminating a costly ground test procedure and 
potentially providing more accurate results than could be obtained from ground based calibration The 
capability to provide ^active sensor calibration” in an operational environment has potential applications, 
beyond the RSRA load cell problem 


1 



NOMENCLATURE 


Symbols 

u . Input 

y .. Output 

... . Number of rotor blades 

0. . . Rotor rate (rad *ec ) 

// . Vector of 6 applied hub loads 

H .Estimates of applied hub loads 

T. ... Vector of 7 loud cell measurements 

J . Vector of 6 inertial loads 
© . Set of unknown parameter* 

0. . -\-pnori estimate of parameters 

© Post-calibration estimate of parameters 

q Set of 1 1 load cell deformat to ; angles 

q Pre-identification estimate of a 

q ... Po-st-identification estimate of a 

it Longitudinal rotor shaft tilt 

5 . Transformation matrix for load ceils 

R . .. Transformation matrix for hub loads 

£*{} Expected value operator 

i7 Vector of measurement noise components 

Q t . Co-variance of measutement noise components 

rj . . Vector of process noise components 

Q, v . Co-variance of process noise components 

x . Vector of states to be estimated 



x . . Estimate of stale vector 
F . Djnamic coupling matrix for states 
. . Transition matrix for states 

n ... . Index of measurement update times 
At Time between measurement updates 

J x . . Partial derivative of inertial loads wrt states 
J „ Partial derivative of inertial loads w r t deformation angles 

P. . Co-variance of state estimation error 
K Kalman gams 

l\ Maxumm Likelihood cost function 

W , . Measurement error weighting matrix 

.V, Number of samples used for identification 

tj. . Error in estimate of inertial loads 

c H Error in estimate of applied hub loads 

jj . Frequencv 

A, ByCy D f £\ F % G Individual load cell outputs (Fig A 1) 

A # , Y % Z f L, Mj S . Force and moment components (Fig A 2) 

6 Position of C G from shaft attach point 

h . Position of hub from C G 

/ . Position of C G from aft load cell attach points 

c . Lateral offset of forward lateral load cell from centerhnc 
/ . . Lateral offset of aft lateral load ceil from centerline 

w Lateral distance between vertical load ceil attach points 

d. Longitudinal distance between vertical load cel! attach points 
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O.0.V 


Load cell deformation angles about x,y and z axes respective!} 



a ... Linear acceleration at C G 
/ . Moment of inertia about C G 

m .... EfTecme mass of rotor/ transmission engine system 
A/ . Mass matrix 

p.< 7 .r . . . Rotational rates about x.y, and z axes 
Q t . Total applied engine and tail rotor drive torques 
T... . Gyroscopic coupling coefhcients 


Subscripts 


A . Accelerometer 

i , j Row and column indices respects el\ for sectors and matrices 

0 Low frequency component 

1 One re\ component 

A . A* rev component 
m Measured data 

j. y,c Component for x y*z axis 

/ Inertial load 

H Hub load 

C ‘Total load 

T Load Cell 

a. 6 c.d,c./,y Attach point of eaui load cell 


Superscripts 

- Pre-measurement estimate 

- Post-measurement estimate 
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PHASE I: DEVELOPMENT AND VALIDATION OF AN ANALYTICAL 
MODEL OF THE ROTOR/TRANSMISSION/LOAD-CELL SYSTEM 

Background 


Previous Methods 

The initial attempt to determine the load cell response to the applied rotor loads involved applying 
static loads at the hub and measuring the resulting load cell response A least squares regression approach 
was used to identif> a coefficient matrix relating the seven load cells to the six hub loads A six element 
bias vector was also estimated It was found 2! that the coefficient matrix \aned as a function of the 
applied multiple axis load This indicates either a nonlinear dependence of the load cell response to the 
applied rotor loads or an unmodeled independent variable in the regression 

The next calibration was a ground based shake test in which a pair ofinertial shakers were mounted on 
the rotor hub to apply dynamic loads at specified amplitudes and frequencies During this experiment, 
the RSRA was suspended from the hub, so the static loading was the same for all tests As a result, the 
nonlinear variation of the relationship with static loads, as observed in the static test, should not have 
been present in the shake test The Force Determination Method 4 was utilized to estimate applied 
rotor loads from a variety of sensors around the aircraft. This method first identifies the transfer functions 
between the applied rotor loads and sensors at various points on the aircraft and then identifies the applied 
rotor load from a least squares fit to the transfer functions and measured sensor responses The results 
were unacceptable because the identified transfer functions varied with the magnitude of the applied load, 
hence they could not be used to estimate the applied load without some e> tensive calibration procedure 
Again, the variation wuth applied rotor load indicates either a nonlinear dependence or an unmodeled load 
source 

Proposed Method 

Since the sensors utilized in the Force Determination Method included numerous accelerometers and 
strain guages mounted on the fuselage, transfer functions of these sensors will be affected by any nonlinear 
dynamic .behavior in the fuselage. This effect complicates the use of sensors on the fuselage to determine 
applied rotor loads The RSRA was designed to use load cells to isolate applied loads from different 
sources, such as the mam rotor, tail rotor, engine and wings The proposed approach takes advantage 
of this concept by treating the rotor transmission as an isolated system with externally applied loads 
from the load cells and the rotor (Fig 1) The applied rotor loads are then measured from a force and 
moment balance using measured load cell loads (7) and inertial loads (J) as derived from transmission 
accleration measurements In order to utilize this approach a model is required that relates applied 
external loads on the rotor/transrnission system to the resulting forces and moments at the center of mass 
of the svstem This model is derived analytical!) from physical principles using the known geometry of 
tiie rotor /transmission svstem Potentially uncertain parameters in this model are explicitly represented 
to provide the capability to calibrate the model 

The advantages of this approach predominant!) arise from the phvsical insight obtained in using an 
analytically derived model With such a model, sensitivity anal) sis and physical judgement can be used 
to select the most appropriate set of available parameters for calibration Consequently fewer parameters 
need be calibrated than when no ph>sical insight used In addition, the parameters to be calibrated 
now have a physical interpretation so the validity of the calibration results ma\ be assessed The model 



should initially be demed to be as simple as possible If it can not adequate!) explain the observed 
experimental behavior with a phvsicallv reasonable set of parameter values, it ma> be expanded to include 
additional effects as required It is important fkat all major effects be identified and incorporated into the 
model before calibration of the parameters is attempted, or the parameter values will compensate for the 
unmodeled effects as best they can and achieve phvsically unrealistic values in the process 

Objective and Approach 

The objective of this pha^e of the ^tudv is to derive a simple dvnamic model of the isolated ro- 
tor/ transmission swem and test its accuracv with experimental data The approach is to first derive 
a simple model of the rotor transmission s)*tem. treating it as a linear, rigid isolated body Known or 
assumed values are u^ed for all parameters of the derived model The model is then applied to test data 
to determine its accuracv *f it appears that calibration can further improve the accuracy, the appropriate 
parameter set will be selected and calibrated If the accuracy appears limited bv an unmodeled effect, the 
model will be expanded to evaluate potential sources of the unmodeled effect 

Rotor/Transmisssion Model 

The arrangement of the load cells below the rotor transmission system is shown in Fig 2 A detailed 
description of this svstem. including all inertia contributions is given in Reference 2 There arc seven load 
cells, four of them mounted vertical!) at the corners of the transmission mounting plate, two mounted 
lateral!) along the fore and aft edges of the mounting plate and one mounted longitudinally at the for- 
ward edge of the mounting plate An inertial load vector, J , is located at the center of gravity of the 
rotor/transmission system and the applied rotor load vector (//) is located at the rotor hub (Fig, 3) The 
rotor hub is located at the end of the rotor shaft which is tilted forward at an angle of 2 degrees The 
load cells are connected to the transmission and the fu-elage bv spherical bearings 

The proposed approach is to estimate the applied hub loads from a force and moment balance of the 
external loads and the inertial loads In order to accomplish this all externall) applied loads must be 
transformed to the center of mass, where the martial load* act Bv treating the rotor transmission as an 
isolated svstem, the load cell force" are considered a measured, exiernailv applied load on the System A 
6x7 matrix (5) is derived that transforms the vector of ^even load ceJT loads (T) at their attach points 
to a set of 6 load components at the center of mass A 6\6 matrix {R) is derived that transforms the 
vector of applied rotor loads (//) at the hub to the center of mass An inertial load vector (J) at the 
center of mass is derived from measured accelerations and assumed inertial parameters. Using the derived 
matrices, a force and moment balance at the center of mass results in a set of six simultaneous equations 
which ma\ be written m matrix form as 


J-ST-RH ~ 0 (11) 

A detailed description of these vectors and matrices are given in Appendix \ 

The assumptions used in deriving the matrices and Eq 1 1 are that the rotor/ transmission system 
is a rigid body and that there is no friction m the load cell bearings These assumptions were made to 
simphfv the initial approach Both non-rigid bod> effects and friction m the load cell bearings could be 
added to the model if that appears to be required in order to match the experimental data 


G 



The transformation matrices 5 and R and the inertia! load vector J were all derived to explicitly 
contain all potentially uncertain parameters of the system so that any subset of parameters can be selected 
for calibration, as required to improve the model accuracy The parameters exphcity defined in the mode! 
are 

a) All distances and angles required to define the resultant moment arm from the center of mass to load 
application points 

b) All angles and magnitudes required to define load ceil load components acting on the rotor/ trans- 
mission system 

c) All mass properties required to determine inertial loads from measured rates and accelerations 

Test Conditions 

Having derived a model, the next step is to evaluate it with experimental data Both static and 
dynamic ground test data are available The dynamic data generated by the shake test were chosen since 
it would provide a more rigorous test of the model structure than the static test data The dynamic 
data are not, however, suitahle for testing the calibration procedure This is because the same static load 
condition exists for all dynamic tests and parameter variations are mostly dependent on variations in the 
static loading Once the model structure has been validated with the dynamic data the model can be 
applied to the static test data to evaluate the calibration procedure 

The test data set selected was a frequency sweep from io to 18 5 Hz in the applied rotor load. This 
frequency range was chosen because it contains the rev frequency, and identification of applied rotor 
loads at this frequency is of special interest Transfer function data were generated from the raw test data 
by a harmonic analyzer for four levels of applied load in each of the three axes Eq I 1 was then used to 
generate the applied load estimate from t he t r ansfer function data Eq I 1 was processed with the values 
of all of the model parameters set to their design specifications 

Error Analysis 


Total Estimation Er^ror , 

In a controlled ground test environment, the actual values of the applied rotor loads are known so 
the total error in the estimate is readily obtained Given the vector of measured load cell readings T m , 
and the vector of inertia! loads derived from accelerometer measurements J n , Eq 11 may be used to 
estimate the applied rotor loads as 


ft = -R- 1 Jm-STl 


(1 2 ) 


The total error in the estimate is obtained by subtracting the known values of applied rotor load from 
the estimate of Eq 12 to get 


e H = 8 - H = -/T 1 - H 

• e» = - Sf m - Ril (1 3) 


t 



The available shake test data was in the form of transfer functions That had been generated from the 
raw data by a harmonic analyze^ In order to utilize this data with the proposed model it was necessary to 
transform the model to the frequency domain and write it in terms of the transfer functions Transforming 
Eq 1 3 to the frequency domain gives. 


?»(-) = RHU) ( 14 ) 

If onI> a y-axis rotor load. is applied, Eq 1 4 maj be written in terms of the transfer functions 

as 

*■#/(-) = -R-' 0 5 ) 

where and are vectors of transfer functions of the inertial and load cell loads with respect to the 

y- axis rotor load Eq 1 5 was used to evaluate the total error in the estimate using the available transfer 
function data for y-axis excitation Corresponding expressions were used to evaluate the estimation error 
for x and z axis excitations. The cfTect of applied moments at the rotor hub u<js not evaluated in this 
stud} 


Numerous potential sources of error are present in this system The evaluation procedure was to 
examine the total error and attempt to categorize it into the potential sources Once identified, the 
sources would be modeled and included in the system Most of the error will probably be atttnbutable to 
one of five sources 

1) Systematic errors in data collection 

2) Unmodeled static and dynamic effects 

3) Unmodeled external loads 

4) Random errors in sensors 

5) Incorrect parameter values 
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Systematic Errors in Data Collection 

One source of error is the use of transfer function data Since this is treated as raw data in this study 
any errors in the identification of the transfer functions would propagate through the proposed approach 
The recorded time domain data should be reprocessed by the harmonic analyzer to provide onl> Fourier 
transformed data, not transfer functions 

Angular accelerometers and rate gyros were not available on the rotor transmission system for the 
shake test Since there is no wav to obtain angular acceleration data, the approach taken is to assume 
its values are negligible and see how good the estimates are with this assumption Some justification 
for this assumption comes from comparing the response of two linear accelerometers that are mounted 
on the transmission with a two foot vertical displacement between them The difference in the y-axis 
components divided by the vertical isplaccmerit should give the roll axis angular acceleration The average 
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value of th*s dem ed roll acceleration over the frequency range wa* found to be 09 rad sec ^.supporting 
the low angular acceleration assumption The derived angular acccderation data wa $ » ot used with the 
model because the errors ir the linear accelerometers are such that the accurac) of the derived angular 
acceleration is G rad f sec sec. The derived values arc therefore in the noise level The anguar acceleration 
affects the translational equations since the linear accelerometer is not mounted at the center of mass and 
will therefore be affected b> angular accelerations The moment equations arc affected since the angular 
inertial loads are dependent on the angular accelrraton* The errors associated with neglecting a 09 
rad *e< sec roll angular acceleration arc 12 lbs m the y force and 27 ft lbs in the rolling moment so the 
assumption appears justified 

Umnodolod Static And Dynamic Effects 

The model derived for thi* study was kept deliberate!) simple to facilitate the analysis It can be 
expanded if necessary to account for the observed error Lnmodeled effects with potential!) significant 
impact on the mode! include friction in the load cell bearings flexibility in the rotor' transmission system 
and nonlinearit) in the dynamic response to applied loads The nature of the error signal should suggest 
which of these effects are present Friction and deadbands will be characteri 2 ed by hysteresis in the 
response This effect is more readily observed in sialic data than in dynamic data Flexibility will show 
up as a resonance at some frequency and will result in a phase and amplitude shift between the input and 
output signals Nonlinear dynamics will be readily detected by a frequency shift between the input and 
output data Static nonlineamtes result in parameter variations and are corrected by calibration rather 
than by expanding the model 

Random Errors In Sensors 

Both load cells and accelerometer* have measurement no^e which produces a lower bound on the 
accuracy of t he applied load estimates It is possible to obtain accuracies in limit but this 

requires the u«e of statistical processing technique* such as Kalman Filtering and Smoothing This effect 
can not be reduced by any modifications to the model The elnct of accelerometer and load cell noi^c 
on the applied rotor load estimate is derived in appendix B and u*ed to generate the errors given in 
Table I This table compares the estimation error* in each axi* for four levels of applied load with the 
predicted standard derivation of the *»rro r based on instrumentation noi**» Th***e numbers were based on 
the assumption of independent random errors for each senior with accuracies of 1 percent of full scale for 
the accelerometers and 0 1 percent of full scale for the load cells The total instrumentation error is shown 
along with the components ddue to load cell and accelerometer noi-r 

Table I 




Load cell full scale ~ 25,000 lb , err ~ l^c 
Vcceleromet^r full scale 7 G's , err — l f 7 
Mass = IKK) |b 
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Umnodelpd External Loads 


The derived model will be in error if any external!) applied load is not included If the levels are low. 
then this can be the most difficult source of error to identify This is because it can take on vtrtuallv any 
characteristics and wilt blend in with other error sources. The onlv possibility for detecting this type of 
error is if it is sufficiently large that it can not be logtcallv explained by any of the other error sources. 
Once the presence of an unmodeled load is suspected, the error data mav assist in isolating its source, but 
a thorough examination of the test condtions is usualU required to r^olve this efTect 

Incorrect Parameter Values 

Once the model structure ha* been validated to the fullest extent possible, the remaining errors should 
be due only to incorrect values for the parameters \t this point calibration may be applied to reduce 
this error source. If calibration is attempted before the model structure is adequately determined , the 
parameters will take on whatever values are required to compensate for the model structure errors This 
amII result m physically unrealistic values for the parameters and eotild, in fact, be a test for whether the 
model structure is accurate 

The conventional approach to calibration ha* been to apply least squares minimization of the error 
with respect to the parameters to be calibrated Since the parameters are now imbedded m a model, the 
least squares minimization of the error must be done subject to the constraint that the model equations 
are satisfied This is refered to as a constrained least squares approach 5^ and the algorithm is derived in 
Appendix C for the constraint of Eq l l 

^mce the available data is in the frequency domain, the calibration must be performed in the fre- 
ouency domain This is actually an advantage since the frequency domain transformation has concentrated 
information for the required frequency ranfre into fewer data points than required for a time domain repre- 
sentation Calibration may then be performed with less data in the frequency domain. The lea<t squares 
minimization may be applied to frequency domain data m the ^amc way as it is applied to time domain 
data C The only modification is that the data is organized with the real and imaginary parts stacked 
end to end rather than using the data m complex form This insures that the identified parameter values 
will not be complex 


Results 

Eq 1 3 was applied to transfer function data for the load cells ami the mam rotor gear box accelerom- 
eters to generate the estimation error of the derived model in t he frequency domain The values used in 
the design specification for the reometry and mass properties were used for all parameters in the model 
Fig's 4 and 5 show the magnitude of the transfer function data for the 7 load tells and 3 accelerometers 
for an SCO lb t/-axi< excitation The location of the load cells mav be seen bv finding the correspondingly 
labeled load cell in Fig 2 The accelerometer transfer functions have been multiplied by the system mass 
to produce inertial loads The strong coupling m the system is apparent from the large vertical load 
cell values for a y-a\w excitation The strong correlation in the r-axis is particularlv apparent from the 
sudden drop in r-axis acceleration at the same frequenev (1GC> Hz) where two vertical load cells ( \ and 
B) suddenly assume equal and oppo c ile values 

Figurts G through 11 show the six axis estimation error*- f or y-axis excitation at 300,400,600 and 800 
pound levels to demonstrate the effect of load level on e>um uiun error. The flight values of the applttd the 
applied rotor load at the .V*/rcv frequency should fall within thi* range figures 12 through 17 «how the 
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six axis estimation errors for r, y, and z axis excitation at the 400 pound level to demonstrate the effect 
of excitation axis on the estimation error The errors plotted m (Figures 6 through 17) were obtained 
from Eq 1 3 using nominal parameter values. Table ! shows the average error over the frequency range 
as compared to the accuracy limit set hy the instrumentation noise With the exception of applied rolling 
moment due to y-axis excitation and x-axis force due to x-axis excitation, the average error shown in 
Table I and the frequency plots comparing the three axes of excitation shown in Fig's 12-17, demonstrate 
that the rigid body model of the isolated rotor transmission system produces applied load estimates with 
error levels comparable to the instrumentation noise level 

The pronounced roll moment error from y-axis excitation (Fig 15) is too great to be explained b\ 
parameter errors, sensor errors or angular accelerations The demonstrated lack of significant angular 
accelerations about the roll axis indicates that it is not due to a non-rigid bodv effect It can also be seen 
that there is no frequenev shift (the hump at 1G6 Hz in the error signal matches the hump in the load 
cell data at that frequenev) so the error is not a nonlinear function of the modeled variables The only 
remaining explanation is that the error is due to an unmodeled extend load with predominant effect in 
the roll axis A potential candidate for the source of this load is the drive tram since it would affect only 
the roll moment Measurements of the shaft torques from the engine and the tad rotor are required to 
verify this and could be used to compensate for this effect. 

The unknown roll torque may be reconstructed from Eq 1 3 to give further insight into its source 
Fig 18 is a plot of the phase of the unmodiled torque relative to the applied excitation as a function of 
frequency for four levels of applied \-a\is load h can be that the phase is sufficiently close to 180 
degrees to indicate that the torque is due to elastic rather than damping effects This is consistent with 
the possibility of a binding in the drive tram rotation 

The pronounced x-axis force error from x-axib excitation (Fig 12) is also too great to be explained 
by parameter errors, senior errors or angular accelerations and no frequenev shift is present between 
input and output, so the tffect is not due to a nonlinearity Again the only remaining explanation is an 
unrnodeled external load A potential candidate for this load source is possible binding of a translational 
coupling on the drive tram that is designed to prevent transmission of x-axis loads from the rotor through 
the drive tratn to the fuselage .Measurements of the tension in the drive tram shaft will be required to 
confirm this and to properly account for this effect in the model 

Out of eighteen estimates (six axes of estimates for each of three axes of applied loads) nominal 
parameter values produced estimation errors within the instrumentation noise level for all but two cases 
Roll moment due to y-a\is excitation and x-axis force due to x-axis excitation had errors too great to 
be explained by parameter uncertainties, so calibration is not needed to improve the accuracy m the 
other Mxteen cases and would not help in the^e two cases The source of these unknown loads must be 
determined and cither be measured or modeled as a function of the measured loads before the model can 
be used to estimate roll moments and longitudinal forces 

Conclusions Tor Phase I 


Concluding Remarks 

\n analytical!) derived model with nominal parameter values has been used to estimate applied rotor 
loads from measured load cell and accelerometer data. This approach has also provided a check on the 
consistency of the measured input output data The pre>ence of an unmodeled roll torque due to y-a\is 
applied loads and an unmodeled x-force due to applied x-axis loads have been detected Reconstruction of 
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the unmodeled loads shows that the> ha\e a phase lag close to 180 degrees from the applied load and are 
therefore predominant!) elastic in nature. Rotational and translational binding m the drive tram appears 
to be the most likely source of these unmodeled loads Shaft torsion and tension measurements will be 
required to confirm this and to incorporate these effects in the model 

Specific Conclusions 

The following conclusions are drawn from th.s stud) 

An anal) lifallv derived linear, rigid body model of the isolated rotor/ transmission system with 
nominal parameter values provides estimates with errors within the level of the sensor noise for all cases 
but roll torque due to applied y-axis load and x-force due to applied x-aus load 

Non-rigid bod> effects or nonlinear behav lor can not explain the large errors present in these two cases 
The error must be due to an external!) applied, unmodeled load The presence of such an unmodeled load 
could also account for the apparent nonlinearities detected b) the previously attempted transfer function 
identification If there are unmodeled load sources the transfer functions will assume whatever values 
are required to reproduce the load cell reading* given the measured load* if the unmodeled loads are 
elastic, and therefore dependent on the applied loads then different tnnsfer functions will be identified 
at different applied loads This is pnciselv the effect observed in identif>ing transfer functions from the 
shake test data 

Reconstruction of the unmodeled load from the error signal indicates it is predominant!) elastic 
in nature If ro, it was probably present in the italic calibration and could account for the apparent 
nonlmeanties detected in the static calibration Applung a regression fit to input/output data in the 
presence of an unknown input will result m the identified coefficients taking on whatever values are 
required to account for the unknown input. \s in the shake rest case, if the unknown loads are elastic, 
they will be dependent on the applied loads and the identifnd coefficients will vary with applied load 
This effect was observed in identifying a model from the static test data 

Calibration of the parameters is unnecessary at this point since the design values appear to give 
satisfactory results when the external loads are all properly included. 

Recommendation For Further Research 


The following recommendations are made 

If the unmodeled loads are e } asnc, rather than damping, as the phase lag appears to suggest, then 
the^e loads were probably pre-°nt in the static test and could account for the apparent nonlinearity in that 
calibration The unalvtica! model should be applied to the existing static calibration data to determine if 
the same unmodeled loads appear to have been present 


The upcoming static calibration should be sufficiently instrumented to detect the suspected sources 
of unmodeled loads from the drive train and the analvtical model should be used to process the data 
with sufficiently fast turn-around that newly identified effects can be investigated before the tests are 
completed 
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The approach should be applied to static test data to evaluate th* calibration technique Nominal 
parameter values have done surprisingly well in allowing accurate hub load estimation for the specific 
static load condition present m the dvnamic data Static test experience suggests that this will not be 
* true under all static load conditions Once the model structure is fully validated it should be applied to 
static tvstf data to determine which parameters to calibrate and to determine the range of variation of the 
parameters through calibration 
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PHASE II: IDENTIFICATION OF ROTOR/TRANSMISSION/LOAD CELL 
MODEL PARAMETERS DIRECTLY FROM FLIGHT TEST DATA 

Background 


Overview 

Phase ! of this study concentrated on the development of an analytical mooel of the rotor* trans- 
mission load cell system and on the validation of that model from ground test data It was determined 
that under some loading conditions, unmodcled external loads are present and these loads must be ac- 
count* d for before the model can be further improved Once this is accomplished, the nett step will be 
to vnhddte the model parameters over a range of operating conditions The available data in Phase I 
tended to confirm that parameter values obtained from the design specifications were adequate but this 
**et of data included only one static load condition and the parameters of the model are more like!} to be 
dependent on the high magnitude static loads than on the lower magnitude dynamic loads It is possible 
that the parameter* will vary with flight condition and must therefore be calibrated against static load 
If the dependency, is a simple one. then this mav be readily accomplished by ground based calibration. 
If the dependency is complicated, however, the required ground testing could be costly and the resulting 
calibration curve may be inaccurate over <ome regions of the operating regime 

Modern parameter identification techniques provide an alternative to laboratory calibration by allow- 
ing direct identification of the parameters from input, output data in the operational environment This 
has the advantage that there is no need to characterize a complicated functional dependency in order to 
determine the parameter values The parameters are identified that provide the best fit of the mode! to 
the input output d,ita over a specified data record If the parameters vary slowly relative to the Tength of 
the data record, the identification mas be more accurate than could be obtained from a calibration curve 
of limited order Tne high cost of generating complicated calibration curve;* from ground testing is also 
eliminated 

Phase II of this study examines the use of modern parameter identification methods in identifying the 
parameters of the R3RA rotor transmission load cell model directly from available flight data In order 
to facilitate the study, the data id was generated by a computer simulation using a model with the same 
structure assumed by the parameter identification method but with different parameter values Since this 
approach only identifies parameters and does not question the model structure it can only be used with 
models whose structure has been adequately validated from ground tests 


Identification and Calibration 

Identification methods can be classified into three categories 

I) Equation Error Method In this method a set of measured inputs, u , and outputs y , nr" assumed 
related by a known mathematical model f(u y) = 0 and the parameters of the model are identified 
by a least squares minimization of the equation error /(«. y) T This technique is commonly referred 
to as regression and requires accurate measurements of both the input and output 



2) Output Error Method In this method, the parameters of the assumed mathematical model are 
identified by a least squares minimization of the output error, determined as the difference between 
the output of the model of the system, and the output of the actual system 8 This approach also 
requires measurements of both inputs and outputs. 

3) Maxtmum Likelihood Method In this method, a Kalman Filter is constructed using a dynamic system 
model and an algebraic sensor model The parameters of both models are then identified by a weighted 
least squares minimization of the output error of the Kalman Filter (innovations sequence) 8 . Note 
that if only the parameters of the sensor model are to be identified the Kalman Filter can be used 
to estimate the inputs to the senior model therefore, only the sensor outputs need be measured 
This feature makes the Maximum Likelihood approach particularly relevant to tensor calibration in 
an operational environment 

Calibration methods can be categorized based on the parameter identification technique used The 
following definitions are proposed 

1) Off-line Calibration . This refers to calibration using either the Equation Error Method or the Output 
Error Method in a test environment where both the sensor input and output are known Once 
the sensor parameters have been identified, they are used with the senior model in the operational 
environment to estimate the sensor input (system states) from measurements of the sensor output 
Since the parameters are not identified from data obtained in the operational environment, this 
method is referred to as “ofT-Iine calibration " This is the moM commonly used approach to sensor 
calibration 

2) On-line Calibration This method uses the Maximum Likelihood Method to estimate simultaneously 
the parameters of the senior model and the input to the sensor (system states) Only the measured 
output of the senior and a model of the d\ namics of the system driving the -senior are needed Since no 
knowledge of the input is required, the parameter identifictton mav be performed using data obtained 
in the actual operational environment hence the term “on-line calibration Note that ~on-!:ne“, as 
used here, does not imply “real-time" Calibration by post-experiment data processing still fits the 
definition of on-line calibration 

OfT-lme calibration is well suited to situations where the parameters are constant or simple functions 
of the environment Under these circumstances, it is the preferable approach because it provides a simple 
means of applying the identified parameters to the sensor model to estimate the states in an operational 
environment When the parameters of the sensor model become complicated functions of the system 
states, however, the off-line calibration procedure becomes more involved It is then necessary to relate 
the parameters to measurable quantities, such as the sensor output, through additional mathematical 
models (usually polynomials) the calibration then involves identifying the coefficients of these additional 
models 

On-hne calibration is ideally suited to handle situations where the sensor mode! parameters are 
complicated functions of time or the system states Since the parameters of the sensor model are direct I \ 
identified in the actual operational environment, there is no need to further characterize the parameters 
as functions of measurable quantities Thus far, on-line calibration lias been used principal!; in support 
of parameter identification of dynamic system models This paper will address on-line calibration as a 
technique for obtaining more accurate calibrations of sensor models with variable parameters 
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Objective And Approach 

The objective of this stud) is to apply the on-line calibration method to the problem of calibrating 
the RSRA load cells The rationale is that b> directly identifying variable parameters of the load cell 
model in the operating environment, the need to model their functional dependence on applied loads is 
eliminated and a more accurate estimate of the applied loads at the hub may be obtained The technique 
was evaluated by using computer simulations of the load cell response to applied loads to generate data 
for the on-line calibration algorithm The emphasis was on calibrating the load cells for high frequency 
applied loads (17 Hz ). because a major objective of the RSRA program is to use the load cells to evaluate 
coupling of rotor loads at the A'w revolution blade passage frequenev (the number of blades, .V* , tunes 
the shaft rotational rate, ) The simulation discussed here is necessarv to refine the on-line method before 
committment to flight tests or other experiments. 


Design Of The On-Line Calibration Algorithm 


Approach 

The design of the on-hne calibration algorithm for the RSRA load cells includes the following steps 

1) Define the load cell model In this step a mathematical model relating applied rotor loads to load 

cell readings and inertial loads for the isolated rotor transmission s>stem is derived The mode! must 
be formulated so that its variation with applied loads can be represented by a set of identifiable 
parameters . , 

2) Define a model of the d\ namics of the applied rotor loads Since knowledge of t lie rotor load dv namics 
is used 2 is an alternative to knowledge of the actual loads m on-line calibration, it is essential that an 
accurate representation of these d\ namics be derived The derived model must have sufficient degrees 
of freedom to allow its transient response to be fit to the experimental data 

3) Design a Kalman Filter to estimate the applied rotor loads The dvnamic model of the rotor loads 
and the modei of the load cell response to rotor loads are combined in a Kalman Filter to estimate 
the net applied loads at the rotor hub The accuracv of the estimate depends on the accuracy of the 
parameters in the load cell model 

4) Design a Maximum Likelihood estimator for the parameter* of the load cell model The Maximum 
Likelihood parameter identification algorithm uses the measurement error generated by the Kalman 
Filter to estimate the parameters of the load cell model Thcr updated parameters are then used in 
the Kalman filter to obtain improved estimates of the applied rotor loads 

The interaction of these four modules with measurements obtained from the aircraft is shown in 
Figure 19 The Kalman Filter consists of the dynamic model and the sensor model linked together b> 

feedback loops with Kalman gams. K The function of the Kalman Filter is to generate an estimate. H 
of the system states, H , using available measurements T and J , to be described below These estimates 

are based on assumed values, 0. for the paramer is of the models, © The Maximum Likelihood Method 
then uses a gradient algorithm on the sequence of Kalrr.an f liter measurement errors , tj , to generate new 
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estimates, 0 of the parameters of the models 8 , For purposes of this stud), the RSRA measurements were 
simulated using computer models that were independent of the models used in the calibration algorithm 
The success of the algorithm wa s then determined by observing the convergence of the models in the 
calibration algorithm to the models in the RSRA measurement simulation 


The Load Cell Model 

The load cell model used for this studv is the model derived in Pha^e ! and described in detail in 
Appendix A The parameter set to be identified vvilf be limited to the load cell deformation angles for 
purposes of this studv, although any of the geometric or mass properties parameters included in the 
derivation in Appendix A could be identified The parameters will be identified from the error between 
the measured inertial loads, as obtained from the accelerometers and mass properties, and the estimated 
inertial loads, as obtained from Eq 1 1 as a function of the parameter set and the applied rotor Ipads 
The model, as given by Eq 1.1, is therefore solved for the inertial loads. J , to give an equation of the 
form 


/= -S(d)f - RH *= .V(0,Q„) (2 1) 

where J is the vector of the six components of inertial loads acting at the center of mass, T is a vector 
consisting of the seven resu .ant load cell loads, and H is a vector of the six components of applied loads 
at the rotor hub. The measurement noise vector, u . is added to characterize the stochastic nature of 
the inertial loads measurement. The noise is assumed to be normal!) distributed with zero mean and 
power spectral density Q v The S matrix is a 6 x 7 geometric matrix that transforms the seven load cell 
resultant loads at the load cell attach points to a ?et of six component loads at the center of mass Note 
that this matrix is a function of the fourteen arbitrary load cell deformation angles dr. that we will use 
as the parameter set to be identified The R matrix is a 6 x 6 geometric matrix that transforms the six 
components of applied loads at the hub to six components of applied loads at the center of mass of the 
rotor/transmission system A detailed description of these vectors and matrices is given in Appendix A 

Equation 2 1 will be used as the sensor model, relating load cell measurements. T to applied rotor 
loads, H Uc have therefore derived a sensor model where the effects of fuselage deformations can be 
represented bv the set of 14 load cell deformation angles Calibration of this model consists of identifying 
these angles Either off-line or on-lire calibration techniques can be used to accomplish this, but since 
the angles are complicated functions of the applied rotor loads, an off-line calibration procedure would 
require modeling and calibrating this functional dependency The on-line calibration procedure allows 
direct identification of these angles in the operating environment in which they will be used to estimate 
the hub loads, consequently, there is no need to model their functional dependenc) on the applied loads 


The assumptions used in the derivation of this model are 

1) The rotor/transmission svstem is a rigid bod), therefore, all changes in the load cell geometry are due 
to relative fuselage deformations Th*s assumption allows the net effect of the fuselage deformation 
on the load cell geometry to be completely characterized bv a set of two orthogonal rotations of the 
load cell axis about each attach point The relative geometry of the applied loads at the hub and 
inertial loads at the center of mass do not change 




2) Applied loads at the load cell attach points are directed along the load cell axis and no moments 
are supported at the spherical bearings This esscntiall) assumes no friction in the bearings and is 
contrary to experimental evidence A significant level of h>stercsis has been observed due to bearing 
fncMon This effect can be incorporated into the model at a later time but will be neglected here for 
purposes of evaluating the on-line calibration procedure 

A third assumption not inherent in the model but u«ed in the on-line calibration algorithm, is that 
the load cell deformation angles vary slow I) so they mav be treated as constant parameters over the data 
sets used for identification 


A Dynamic Model of the Applied Rotor Loads 

A key element in on-line calibration is the use of a dvnamic model of the system driving the sensors, 

It is this model, in connection with the time varying data record that allows calibration with specific 
knowledge of the sensor tnput s The application addressed here requires the modeling of the d>namics 
of the applied rotor loads stnee this is the input to the load cells Modeling the rotor load dynamics J 

across all frequencies would be a formidable task, requiring accurate knowledge of the dynamic load j 

response to pilot inputs and changing flight environment There are frequencies, however, at which the 
dynamics of the rotor loads are well known. In a rotorcraft. strong vibratory loads are present at the one/ 
revolution frequency and at the .V s 'rev frequency and its harmonics, wdiere .V* is the number of blades \ 

of the rotor (currently AY - 5 for the RSRA) These vibratory loads may be modeled using second order 
systems with zero damping and resonant frequencies at the prescribed values By filtering the load cel! 
and accelerometer data to pass only these frequencies, the d>namic content m the filtered data will be j 

accurate!) represented by the specified models This dynamic content is essential to the identification of j 

even quasistatic parameters since there are more parameters than equations and the problem is therefore t 

statically indeterminate [' 

»• 
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By filtering the data to pass only selected frequencies with v ibratory loads, the pilot input is eliminated * 

from the data Pilot input is normally used in parameter identification to insure that the parameters are I 

identifiable by providing sufficient degrees of freedom in the dynamic response and by exciting all modes j 

of the dynamic system Since the parameters to be identified are all in the algebraic sensor model there I 

arc no specific modes to excite, but the input to the load cells must have at least as many degrees of * 

freedom as there are parameters to be identified In this case there are fourteen unknown parameters If i 

we model the vibration at one frequency as a second order svstem for each component of hub load, the j 

dynamic response will contain 12 independent components (a sine and a cosine component for each of the 
six hub loads) We therefore need to pass two frequencies with vibratory loads so that the dynamics will * 

contain 24 degrees of freedom and insure identifiabihty for the 14 parameters Using this approach makes 
the parameter identification process independent of the pilot input allowing the calibration to proceed in | 

parallel with any chosen flight maneuver j 


Since two v ibration frequencies must be passed for adequate sensor excitation, the two lowest frequency v , 
vibration modes associated with the rotor wall be used These are the one rev and the AY/ rev frequencies i 

The dvnamic models for the one rev and the AY, rev vibrations are given by the following equations 
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where 0 is the rotational frequency of the rotor and r)i and rj* are normally distributed'random variables 
with zero mean and power spectral densities Q ni and Q ns respectively Both H\ and //? are six component 
hub load vectors and the total hub load vector present in the filtered data is given by 

H = H t - /7a (2 4) 

While Equations 2 2 through 2 4 are not physically based models of the vibratory loads, they provide 
a mathematical representation that can be used to fit the dynamic content of the filtered data By 
specifically. -<&aracteriz»ng the frequency characteristics, a more accurate fit to the data can be obtained 
for a model of given order 


Equations 2 2 and 2 3 may be written m state variable form as 



■W, ‘ 


r o 

/ 

0 

°i 

■ Hx ■ 


- 0 - 

d,dt 

H , 

_ 

-n : / 

0 

0 

0 1 

H , 


'U 



0 

0 

0 

/ > 

//% 


0 


■ Hs. 


. 0 

0 


oj 

JU. 


.n \ - 


Defining a state vector, x, by 


we may write Equation 2 5 as 


f = [//. 
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The \alue of z at time n x(n) mat then be computed from f(n - 1) b> 


(2 5) 


(2 6 ) 


f(n) = <f>(Ju)i(n - 


At 

1) - J 4>(Af - r)rj(r)dr 
0 


(2 *) 


where the time inter' a) from n - l to n is Af and 

$>(Af) = e FM 


(2 9) 


Design of the Kalman Filter 

Having defined a dynamic model of the rotor load^ and a geometric model of the response of the load 
cells to applied rotor loads and inertial loads a Kalman Filter to estimate the rotor loads may be derived 
The sensor model Equation 2 1 will be used as the measurement equation but the load cell readings, T. 
cannot be used au the measurements This is because Kalman filter must estimate the measurements from 
its state estimates and use the error between the estimated and actual measurements to update the state 
estimates There are seven load cells but only six mea- surernent equations m Equation 2 1 As a result 
a unique set of estimated load cell readings, 7\ cannot be generated from the six state estimates using the 
six measurement equations The solution to this problem is to use the inertial loads because they may 
be uniquely estimated from the six estimated rotor loads and the :*e\en known load cell loads using the 
six measurement equations of Equation 2 1 It should be noted that the inertial loads cannot actually be 
directiv measured but must be computed from available accelerometer measurements and assumed mas:* 
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properties The accurac) of the assumed mass properties and accelerometer responses could significant!) 
affect the results 

Equations 2 !, 2.4, and 2 6 ma> be combined and written in the form 

/= -S(d)f - R H x - = J(a.f,x) (2 10) 

The Kalrnan Filter 5 ma> now be written using Equation 2 10 as the sensor model and Equation 
2 7 as the dynamic model Between measurement updates the state estimates, x, are propagated by the 
deterministic part of Equation 2 8 as 

£~(n) = - 1) (2 11) 


where the superscripts - and — refer to estimates before and after the measurement update at the specified 
time 


The post-measurement estimates ma> then be obtained from apriori estimates at the measurement 
update time, n, by 

x (n) = x (n) - K(n) J m (n) - J(a.f{n) x (n;)! (2 12) 


where J m is a vector of 6 inertial loads derived from measured accelerations Note that the stochastic 
part of the dynamic model. Equation 2 8, is simulated by the measurement error feedback in the Kalman 
Filter, Equation 2 12 The Kalman gains /t(n). are given by* 

A'(n) = P~(n)J z T J 2 PJj — Q v \ (2 13) 

where 

J t = - R 0 R 0 (2 14) 

p-(n) = 4>P~(n- 1)<P T -Q„ _ (2 13) 

and 


P~(rt) = ,1 - K(n)J x P~(n) (2 16) 

P“ and P^ are the apriori and post measurement variables of the state estimation error, respective!) 

Since J z is constant and K is independent of the identified parameters, q, Equations (2 13) through 
(2 16) will reach a stead) state condition that is independent of the identification process The algorithms 
can be greatlv simplified bv using the resulting stead) state value of the Kalman gains, K, m Equation 
2 12, thereby eliminating the need to use Equations 2 13 through 2 16 The use of a stead) state Kalman 
filter where the gain is independent of the parameters to *>e identified insures a globally stable Maximum 
Likelihood Identification Algorithm 9* 
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The Maximum Likelihood Estimation Algorithm 

The accuracy with which the Kalman Filler can estimate the applied hub loads depends on the 
accuracy of the estimate of the load cell deformation angles q In the proposed approach these values are 
updated on-line using Maximum Likelihood parameter estimation The Maximum Likelihood estimation 
algorithms are derived by minimizing a weighted sum of the squares of the measurement error of the 
Kalman Filter, (the innovations sequence) with respect to the parameters over a data set of X, samples 
The cost function is defined as 

iV. 

(2 17) 
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(2 18) 


and 


W = l/N,^ej(n)e T j(n) 


(2 19) 


Denoting the apnon estimate of the vector a b> d, and the post-measurement estimate by q, the first 
partial of with respect to the deformation angles is expanded about its value at d to give 


dv(d) av{&) _,d*v(*) 

Setting the desired value at d to zero corresponding to an extremum, gives 
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Taking the first partial of V with respect to o gives 
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where 
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(2 23) 


The dvnamic behavior of the partial derivatives of the aprion state estimates, x . with respect to the 
angles, a is obtained from Equations (2 11) and (2 12) as 
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The partials of the estimated inertia) loads with respect to the angles J ai are given by 
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The second partial, to first order, is then 
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The parameter estimates are obtained from Equations 2 21 through 2 28 using the measurement error 
sequence, tj , generated bv the Kalman Filter 


Operational Approach 


The data is first filtered to pass onlv the two specified frequencies The Kalman Filter is then applied 
to the data and allowed to reach steady slate before the Maximum Likelihood algorithm begins operating 
on the measurement errors generated b\ the Kalman Filter The reason for waiting is that the error m 
the initial rotor load estimate will bias the estimates of the load cell angles By waiting until the transient 
due to this initial error decays we can eliminate this source of bias without having to add the initial 
conditions to our set of parameters to be identified This delav is necessar) only on the first data used for 
identification If subsequent data sets are contiguous, the initial estimate of the hub load can be set to the 
final estimate of the hub load on the previous data set and the initial error will be small The statistics of 
the process and measurement noise. Q n and are chosen to generate stead} state Kalman gains from 
Equations 2 13 through 2 16 that provides a high filter bandwidth (short time constant) compared to the 
frequencies of the vibrations being tracked 1 1 0 

The Kalman Filter designed for use in calibrating the load cells will estimate onh the one/ rev and 
•W/rcv components of the hub loads If t u c low frequenc> response is also desired, it is necessary to pass 
the unfiltcred data through a Kalman Filter that utilizes the updated parameters An additional 6 degrees 
of freedom dynamic model must be included in this Kalman Filter to track the low frequenev components 
of the hub loads This additional states can be modeled by the Weiner process 


(2 29) 


Ho = no 

and the total hub loads are then given b\ 


rfo = A (O.Qrj.,) 


// = Hu- H x - //v 


(2 no) 

Simulation Results 

The algorithms wpre implemented as shown in Figure 19, and a computer simulation of the applied 
hub loads and the measured load cell respond wa* u*ed to generate simulated measurement data for the 
HSRA Again the problem of having to generate seven unique load cel! readings uMng the specified hub 
loads and inertial loads with the six measurement equations was encountered For purpose* of simulation, 
tlie problem wa* solved bv adding a seventh constraint equation whirh e*.seiuially minimized the strain 
tnergv in the hud ct IN Thi> approach yielded reasonable Mmulation values ofload cell readings for the 
applied hub loads The physical basts for this constraint equation was too questionable to allow it to 
be utilized in the identification algorithm, so it was only used to generate simulated data with which to 
validate the algorithm. The deform Uion angles in the load cell simulation were set to arbitrary constant 
values and a different -et of values were used with the Kalman Filter to initiate the identification process 
Values of zero degrees were u*ed for the simulation and four degrees for the Kalman Filter 

No inertial loads were simulated for the test since this represents a worst case for the identification 
process Any input to the system from the inertial loads would have added additional degries of freedom 
to the input to the load cells and thereby improved the idcntifiabihty of the load cell angles Note that the 
presence of such additional measurement inputs would reduce the required number of degrees of freedom 
in the vibration model 

As explained in the previous section the initial values of the hub loads and the load cell deformation 
ancles wore set differently in the Mmulation model and in the calibrrtion algorithm to evaluate the abilitv 
of the algorithm to converge to the value* used in the Mmulation Transient's re<uhmg from the initial 
state error in the Kalman Filter estimate will bias the parameter estimates In order to avoid this problem, 
me filter v\a* allowed to reach a steady state condition before the Maximum Likelihood algorithm was 
initiated 

The tune constants chosen for the Kalman Filter were all at 0 1 seconds so a delay of about 0 5 seconds 
was necessary to insure that the response had readied a steady state condition Two frequencies, one rev 
and A\ / rev , were included in the simulated hub load dynamics driving the load cell simulation. This was 
to insure that the load cell input had mfficient degrees of freedom to permit identification of the 14 load 
cell angles It also effectively simulated filtering all but those frequencies in real test data A sampling rate 
of CG samples, second was used to insure adequate reconstruction of the 17 Hz. 4 V% revolution frequency. 

Figure 20 is a plot of the error m the estimate of the longitudinal deformation angle of one or the 
vertical load cells \ delay of -eventv samples m implementing the identification was used to allow the 
filler to reach slcadv state A sequence of 10 samples was then u*ed to perform the estimation Due to the 
nonlinear nature of the measurement equation, multiple iteration^ over a data ** t vere required to obtain 
convergence for the parameter et-turate^ On each iteration the sevemv sample delav was repeated and 
the last 10 sample* were used to estimate the load cell angle* Excellent convergence was obtained within 
four iterations over the same data set Figure 21 is a similar plot for a lateral load cell deformatiori angle 

Figure 22 compares the actual y-a\is hub load with estimates obtained using the initial and final 
parameter estimates The e=tirnat^> were obtained from a Kalman Kilter with a Weiner process model of 
the additional low frequency states The Kalman Filter um d a simulated unfiltercU data set that contained 



a low frequency component in addition to the two vibrator) components The hub load estimation error 
is significant for the initial parameter estimates but is considerably improved when the final parameter 
estimates are used A similar result is seen in Figure 23 for the applied roll moment at the hub The 
computer simulation demonstrates dearlv the potential for performing on-line calibration of the load cells, 
using data obtained in the flight environment to simultaneous!) calibrate .the load cells and estimate the 
hub loads Figures 22 and 23 demonstrate the sensimit) of the load estimation to errors ta the load cell 
angles The load estimation error is significant!) reduced when acurMe estimates of the parameters are 
used 


Identification From Ambient Vibrations 

The demonstrated capability to identify parameters from ambient vibrations j* significant in that it 
allows identification to proceed in an operational environment without the need to superimpose special!) 
designed input signals to facilitate the identification The method ma\ be applied to other problems 
beyond that discussed here For any such application the conditions for implementing this approach are 
as follows. 

fj All parameters to be identified must be included in an algebraic representation of the system since 
dynamic models would require excitation ai specific frequencies This does not present a problem since 
any dynamic model can be treated as a set of algebraic equations if all stales and their derivatives 
are explicitly measured or estimated at each point in time 

2) Ambient vibrations must include as manv degrees of freedom as there are parameters to be identified 
and the dynamic model of the vibrations must also include these degrees of freedom This can be 
accomplished by modeling only the required number of degrees of freedom and filtering the data to 
pass onlv the modeled frequencies 

Potential applications of this approach include the tracking of time varying stability derivatives of 
a flight vehicle during maneuvers that were designed to investigate performance This wcu*d provide a 
direct correlation between design parameters and observed handling qualities A real-time implementation 
of this identification technique would also facilitate adaptive control applications. 

Conclusions For Phase II 


Concluding Remarks 

The purpose of Phase II of this stud) was to analytically demonstrate the feasibility of an on-line 
calibration procedure Rather than use experimental data, the investigation was conducted using dati) 
generated by a computer simulation of the load cells to allow a companion cf the identified model with 
the known simulation model The structure of the mode* to he identified, as implemented in the Kalman 
Filter, was the same as that used to generate the simulated load cell data, but the states and parameters of 
the two models were initialized differently The unknown model parameters (load cel! deformation angles) 
were then identified succe*:>fullv by applying the algorithm to the simulated load cell data The resulting 
estimate of the states (applied rotor loads) also converged to tile values generated by the simulation. The 
model structure assumed for this investigation has not been validated with experimental data and the 
effect of using a mode! structure in the Kalman Filter that does not match the simulation mode! has net 
been investigated 



Specific Conclusions 


The following conclusions can be drawn 

1) For the assumed model, the dynamic rotor-load estimates were found to be very sensitive to static 
load cell geometry deformations The rotor-load estimates were greatly improved by using identified 
geometric parameters 

2) Calibration of all model parameters was successfully accomplished using only measured load cell and 
inertial loads data 

3) The identification of the model parameter* was accomplished using only rotor vibrations at the two 
lowest frequencies (one rev and A* /rev) as the source of excita- non. No specially designed external 
inputs were required 

4) The identification was accomplished with the inertial load data set to zero This represents a worst 
case situation since any independent degrees of freedom in the inertial loads would provide excitation 
that would enhance the identification 

Recommendation For Further Research 

The following simulation-ba^ed investigations will be necessary to fully qualify the on-line calibration 
procedure for application to flight data 

1) The rapabibtv to track time-varving deformation angles mu*t be tested using simulated load cell data 

1!) The effect of senior nui«« and model structure error on the identification must be determined 

3) The effect* on performance of varying the modeled vibration frequencies, the sampling rates and the 
Kalman T liter time constants should be investigated 

A practical implementation of this mcthodologv will require a validated mathematical model of the 
load cell response to applied rotor loads The accuracy of this technique is a No dependent on the accuracies 
of the assumed mass properties and accelerometer responses. All of this must be determined from suitable 
ground tests before the on-line calibration approach can be applied to flight data 
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APPENDIX A: LOAD CELL MODEL DERIVATION 


A diagram of the undeformed load cell geometry is shown in Figure A 1. The load cells are mounted 
to the fuselage and to the transmission base by spherical bearings It is assumed that the transmission 
base does not warp so all changes in the load cell geometry are due to deformations in the load cells 
themselves or in the fuselage mounting points With this assumption it is possible to completely model 
the load cell response using a general three component representation of the reaction force at each attach 
point on the transmission base, no knowledge of the fuselage deformation is required 


Three plane views of the load cell geometrv with the three component reaction force representation 
are shown in Figure A 2 Inertial loads including gravity and the drive *haft torque, are assumed to be 
concentrated at the transmission center of gravity The rotor loads are concentrated the hub Taking 
the sum of the forces and moments about the cei-ter of gravity gives 


A\ = X t - Xu -G x ~ E x ~ F x ~ A x - B x -C x - D x = 0 


K - Yi - Yu - C y - Ey - F y - A y - By - C„ - = 0 


(.4 1) 


~ Z / — Z yi G t — E z F z — A* B £ — C t D £ — 0 
L c ~ Li - Lu - Yiihz - (A, - D t )u. 2 - E : t - {B : - C 2 )w 2 - F z f 
~( By -C y - F y -Gy - E y - A y - D y )b z = 0 

A/c = A t t - Mu - X H hz - Z !f hx - (B, - A m -E z - G x )(d - C) - (C T ~D X -F X ~B X - A x -E x - G x )bz 

-{C x + Dm- F t )C=0 

X c = Xi - Xu - ~E v ~Gy~ By)[d - L) - (C y - F y - D y )t 
~{A Z - D x )w/2 - E x e - (C r - B : )w/2 - F x f = 0 

Using the transformations. 


A x = Acos(o a )sin(0 a ) 
A y = -Astn(o a )cos(0 a ) 
A t = Acos(o a )cos(0 a ) 
B x = Bcos{6b)sin[9 b ) 
By = -Bstn(o&) 

B z = Bcos(o>.)cos(^<.) 
C r - Ccos(o jMn(0 c ) 
C v - -CMM(Or) 

C, = Ccos(o.)co>(^ c ) 

D x = Dco«»(0j)*in(0j) 


(A 2) 


D v ~ - Dsm{os) 


D t - Dcos(o i)cos(9 i) 






u 



E x = -£cos(o < )5in(y f ) 
E y = Ecosfojco^v,) 
E t = E5in(o f ) 

F x = -f'co5(o/)sin(v/) 
F y = Fcos(o/)coj(t*./) 
F, = F$tn(o/) 

G z = Gco$(0 ; )co*(v,) 
G y =: Gcoa'(0 y )sin(i'j) 
G, -= -Gsinffl,) 

the equations become 


J - RH - ST = 0 


(.4 3) 


where 


/*■ = X h Y h Z h L,.M$.Xi 


(A A) 


ind 


r-v/i 


-rn T <l z 

}•/ 


-m y a i -T y q 

2/ 


-m,a. 

L, 


- IzP ~ <?! ~ Ti? 

V; 


-/»?- T \ t q 

L.VfJ 


-l,r 


(.4 5) 


la equation A 3, a It oj/ . and a 9 are linear accelerations in each direction measured at the rotor' 
transmission sy stem center of gravity and p q and r are rotational accelerations Because of the nonrigid 
engine mountings the engine contributions to inertial forces are not equal in all directions when measured 
at the system center of gravity This effect can be adequate!) modeled by assigning different values to the 
total effective mass of the combined rotor transmission engine s>stem in each direction. rn z ,m y , and 
m t There are also a few minor error terms not given here that are discussed full) in Reference 4 / x , / y . 
and /, are moments of inertia, related terms in the cross-products of rotational rates are neglible Q t is 
total applied engine and tail rotor shaft torque Gyroscopic coupling forces due to engine and transmission 
rotational moments of momentum are represented by the coefficients f. with subscripts for the appropriate 
a\is 

The hub forces (//) and load cell readings (T) in Eq A 3 are given by 


H- A'// A///. V w ^ (A 6) 

f = A.BX.D E.F.G (A 7) 

The geometric transformations from the applied loads to the center of mass (R and S) are given by 


R = 


r i 

o 

o 

o 

-hz 

0 


0 

1 

0 

hz 

0 

0 


0 

0 

I 

0 

hz 

0 


0 

0 

0 

0 

1 

0 


(.1 8 ) 


and equation (A 0) (on the next page) in which C and S denote cosine and sine functions, respectively 


h 
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APPENDIX B: THE EFFECT OF SENSOR NOISE ON LOAD ESTIMATION 
The analytical model of the rotor 'transmission system derived m the text has the form 


/ ~ ST - RH = 0 ' . (B.l) 

Given the measured load cell readings, 7 m , and inertial loads derived from accelerometer measure- 
ments J m , the applied rotor loads, //, are estimated from Eq, B 1 as 

‘h = -/?-> (D 2) 

The measured load cell and accelerometer values mav be written in terms of their true values and a 
random measurement nui>e component as follows 

“m = a — I'A. »A - X(Q-Qa) ' (S3) 

T m = f~ »T - A’(0 ,Q T ) ( B 4) 

Writing the derived inertial load vector.J m , as the product of an inertia matrix, A f, and the ac- 
celerometer measurements ci m , gives 

Jm = A /«„ = Ua-\1» A = /- A \v A (B 5) 

where J is the actual inertial load vector 

Substituting Eq's B 4 and B 5 into Eq B 2 gives 

H = -R-'[J-Sf:-u H (B6) 

where 

u j{ — — A — S * (B 7) 

Substituting Eq B 1 into Eq B 6 then gives 

fi = H~P n ^ = .V(OQ w ) (B S) 

From Eq B S it is seen that as given b> Eq (B 7) is the combined effect of the instrumentation 

errors on the applied rotor load estimate,// This error represents a lower bound on the accuracv of the 
estimation that is attainable without applying statistical processing, such as Kalman filtering or smoothing 
The covariance of this error. Qh , may be computed from Eq (B 7) using the known covariance of the 
instrumentation errors. Q * and Qr , sis follows 

Qn = E{uh^h) =£{ -/r 1 A/iM - Si ~T. 

-R-' \te*-Sffr T ) (B 9) 

Qn =!r'ME{vrf}M T {R-') T - S E{u T ^}M T R~' T 

- R- l ME{P^)S T {R- l Y - R-'SE{v T i7$}S'R- 1 (B 10) 

The sensor noise components are assumed independent so the term E{i/ a I't t } is zero Substituting the 
known covariance matrices of the sensors for the other expected value terms gives 

Q„ = /r 1 \1Q a U t - SQ t S t (R-') t (B 11) 
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APPENDIX C: CALIBRATION ALGORITHM 


The applied rotor load estimation error is given bv Eq 1 3 in the text as 

e„= -ir'[J n -Sf m - RH (Cl) 

The coefficient matrices, R and S contain geometric parameters of the model and the derived inertial 
load-, vector, J m . includes inertial parameters of the model A subset. ©, of the parameters is selected for 
calibration and the error is treated as a function of tho^e parameters \ cost function is written in the 
fornr 


r( 0 ) = 1 2y<„|e.n) r ll -‘e«(e n ) (C 2) 

n - 1 

where W is a weighting matrix given bv 


IV = 1/ iV, ^ e//(0, n)e7/(©, n)i T (C3) 

n= 1 

Denoting the aprion estimate of the vector © by © and the postmeasurement estimate by ©, the first 
partial of V with respect to the parameter 6 is expanded about its value at the aprion parameter estimate 
© to give 


me) = av(e) _ • _ 6 a-v(0) 
a© de a-e 


Setting the desired value at 0 to zero, corresponding to an extremum, gives 


(C- 1) 


c>-V(0) . , cH’(9) r 
de 2 ' dQ 


Taking the first partial of V with respect to © at © gives 


dV{Q) 

dO 


A. 

] Te7/(©.n) r » 

ti-1 


-, dt H (© n) 

de 


The second partial to first order is then 


ga (Q) _ 3eW(0 n) T|v ., 3g- w (8.n) 

de 2 fr, a© a© 


(C.5) 


(C6) 


(C7) 


Once the parameter set has been chosen, the partial derivative of the errors, t with respect to the 
parameters must he generated analytically from Eq (C 1) The parameter estimates are then obtained 
from Eq > (C =>) tu (C 7) using the measurement errors sequence, f a generated by Eq (C l) 










Figure 4 Magnitude of load cell transfer functions for 800 lb applied lateral force 



Figure 5 Magnitude of inertial load transfer functions for 800 lb applied lateral fon 
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MAGNITUDE OF THE ERROR IN Y AXIS, lb £ MAGNITUDE OF ERROR IN X-AXIS, lb 



agnitude of r-axis applied force estimation error for applied lateral forces 



Figure 7 Magnitude of y-axis applied force estimation error for applied lateral force- 




Figure 8 Magnitude of z- axis applied force estimation error for applied lateral forces 
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Figure 9 Magnitude of x-a\is applied momert estimation error for applied lateral force 
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Figure 10 Magnitude of i\\s applied moment estimation error for applied lateral forces 



Figure 11 Magnitude of r-axis force estimation for applied r y and r axis Mid 
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Figure 19 Schematic diagram of on-line calibration 



NUMBER OF ITERATIONS ON 10 SAMPLE DATA SET 


Figure 20 Identification of vertical load cell angle 
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NUMBER OF ITERATIONS ON 10 SAMPLE DATA SET 


Figure 21 Identification of lateral load cell angle 
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